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Noisy saltatory spike propagation along myelinated axons is studied within a stochastic Hodgkin- 
Huxley model. The intrinsic noise (whose strength is inverse proportional to the nodal membrane 
size) arising from fluctuations of the number of open ion channels influences the dynamics of the 
membrane potential in a node of Ranvier where the sodium ion channels are predominantly localized. 
The nodes of Ranvier are linearly coupled. As the measure for the signal propagation reliability we 
focus on the ratio between the number of initiated spikes and the transmitted spikes. This work 
supplements our earlier study [A. Ochab-Marcinek, G. Schmid, I. Goychuk and P. Hanggi, Phys. 
Rev E 79, 011904 (2009)] towards stronger channel noise intensity and supra-threshold coupling. 
For strong supra-threshold coupling the transmission reliability decreases with increasing channel 
noise level until the causal relationship is completely lost and a breakdown of the spike propagation 
due to the intrinsic noise is observed. 

PACS numbers: 

I. INTRODUCTION 

As the fast propagation of action potentials along the axon is of fundamental importance in the nervous system, 
e.g. for the successful evolution to large body sizes of organisms, or the information processing in the brain, the study 
of the propagation mechanisms is of great interest for neuroscientists, physiologists and physicists [![. Since the 
empirical modeling proposed by Hodgkin and Huxley in 1952 0], the neuronal firing dynamics with respect to spike 
generation and signal propagation are successfully studied within this deterministic modeling. In recent days, issues 
relating to the constructive role of noise on these dynamics were addressed 0, H| . The noise-assisted enhancement 
in weak signals transmission, transduction or detection has been investigated, e.g. in the context of noise supported 
wave propagation in sub-excitable media Q or in excitable systems The conductance fluctuation of the neuronal 
membranes which arises from random channel opening and closing can, a priori not be neglected @. Indeed, it 
was shown, that this intrinsic channel noise Q can lead to generation of so-called spontaneous action potentials, to 
synchronization phenomena like stochastic resonance @~tH| and coherence resonance (l2l . [l~3j and to synchronization 
of ion channel clusters [l4[ . 

Even the saltatory spike propagation which results from a highly non-uniform distribution of the ion channel, can 
take benefit of the intrinsic channel noise as we have shown recently in Ref. [l5[. The saltatory spike propagation 
occurs in myelinated axons where the activating sodium ion channels are concentrated at the nodes of Ranvier, which 
are separated by segments sheathed with myelin. This results in a much faster propagation speed in myelinated axons 
than that in un- myelinated axons with constant ion channel density fl6l . [ItJ . 

With this work, we extend our prior study [lj| on the effect of channel noise on the propagation of action potentials 
along myelinated axons. In terms of transmission reliability we are discussing the influence of the coupling strength 
between neighboring nodes of Ranvier and that of strong intrinsic noise. 

II. MODEL 

In order to model the signal transmission along myelinated axons, we consider a compartmental stochastic Hodgkin- 
Huxley model [L5J . Accordingly, each node of Ranvier is modeled by a stochasticgeneralization of the Hodgkin-Huxley 
model, which extends the applicability of the original Hodgkin-Huxley model [2| towards stochastic dynamics of the 
membrane potential of finite-size ion channel clusters 0, [l8[ • Each node of Ranvier couples linearly to its nearest 
neighbors. Consequently, the membrane dynamics Vi at the ith node of Ranvier reads (with i = 0,1,2, N — 1, 
where N corresponds to the total number of axonal nodes of Ranvier): 

C^rVi = 4io„ic(^) + Winter (*) + 4cxt W , for i = 0, 1, 2, N —1, (la) 
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FIG. 1: Sketch of the myelinated axon: Each node of Ranvier is treated within a stochastic generalization of the Hodgkin-Huxley 
model and is bi-linearly coupled to the nearest neighboring nodes of Ranvier. 



with the ionic membrane current (per unit area) within the ith node of Ranvier given by the Hodgkin-Huxley model 
1 

hionUVi) = -G K (ni) {V, - E K ) - G Na (mi,hi) (V$ - £ Na ) - G h (V t - E h ) , (lb) 
the inter-nodal currents 

U{V i+1 -Vi) for i = 0, 

ivnter(i) = < K (Vi-l -V t ) for i = N — 1 , (lc) 

I k {Vi-i — 2Vi + Vi+i) elsewhere 

and the external current stimuli Ii yex t(t) at the ith node of Ranvier. In Eq. ([Taj) , C denotes the capacity of the 
axonal membrane per unit area and is given by C = 1/iF/cm 2 . The coupling strength between next- neighboring 
nodes is characterized by k, cf. Eq. (fTc|). According to Hodgkin-Huxley model the reversal potentials read for the 
sodium current: £jMa = 50 mV, for the potassium current: Ek = — 77 mV and for the leakage: El = — 54.4mV. The 
conductances per unit area are given by: 

G K («i) = 5K ax nf, G Na (m 4 , h) = ™\h, . (2) 

and the constant leakage conductance Gl = 0.3mS/cm 2 . In Eq. ([2j, and g^ x denote the maximum potassium 
and sodium conductances per unit area, when all ion channels within the corresponding node are open. For simplicity, 
we assume that every axonal node has the same kinetics, i.e. the same number of sodium and potassium ion channels. 
So the maximum potassium and sodium conductances g™* = 36mS/cm 2 and g™a X = 120mS/cm 2 are identical 
constants for every node of Ranvier. 

The gating variables nj, m, and hi in Eqs. (|lb|) and ([2]), describe the open probabilities of the ion channel gates 
in the ith node, and undergo a stochastic process which stems from a birth-and-dcath-like process of the gating 
dynamics. The dynamics of the gating variables are voltage dependent, and are governed by the set of (Ito)-stochastic 
differential equations [l8l - [2Cl | , 



A Xi = a x (Vi) (1 - Xi) - f3 x (Vi) x t + ii tX (t) , (3) 

with x = m,h^n. Here, £,i tX {t) are Gaussian white noise with vanishing mean and vanishing cross-correlations. For 
the same node with the nodal membrane size A the non-vanishing noise correlations take the following form: 



<6,m(t) km(t')) = [a m {Vi) (1 - mi) + p m (Vi) mi] 5{t - t') , (4a) 

(ti,h(t) Zi, h (t)) = Wh(V t ) (1 - hi) + p h (Vi) hi] 6(t - t>) , (4b) 

(&n(t)ti, n ( t ')) = ^— [a n {Vi){l-ni) + p n {Vi)ni]5{t-t'), (4c) 

where the ion channel densities are pNa = 60 ^m -2 and pk = 18pm~ 2 . The noise strength is decided by the nodal 
membrane size A which is the same for all nodes. In Eq. [3l the dynamics of the opening and closing rates a x (V) 
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FIG. 2: Deterministic transmission reliability: The dependence of the transmission reliability 1Z on the inter-nodal coupling 
strength k is depicted for the deterministic case. For sub-threshold k < 0.067 mS/cm 2 no spike propagation is observed, i.e. 
1Z = 0. Perfect spike transmission, i.e. 72. = 1, is found for the supra-threshold case k > 0.136 mS/cm 2 . In the intermediate 
range rational numbers k : I for the transmission reliability are found. 



and (3 X (V) [x = m,h,n) are taken at T = 6.3 °C. They depend on the local membrane potential V and read (with 
numbers given in units of [ml/]) [3, [2l|: 

0.1(V + 40) 

amiY) = l-exp{-(T/ + 40)/10} ' (5a) 

/3 m (^)=4exp{-(y + 65)/18}, (5b) 

a h (V) =0.07 exp {- (V + 65)/20} , (5c) 

& (F) = TT r l xqn/inl . ( 5d ) 

1 + exp {- [V + 35)/10} 

0.01Q/ + 55) 

= l-exp{-(T/ + 55)/10} ' (5e) 

(3 n (V) = 0.125 exp {- (V + 65)/80} . (5f) 



III. SPIKE TRANSMISSION 



In order to analyze the transmission reliability we exemplary consider a chain consisting of ten nodes of Ranvier, 
i.e. N = 10 and numerically simulate Eq. (JTJ). By applying a constant current stimulus on the first node only, i.e. 
we set /o.ext = 12 /^A/cm 2 and Ji i0X t = for i = 1, N — 1, we find that action potentials are periodically produced 
that propagate along the transmission line. We define the transmission reliability coefficient 1Z in steady state by the 
ratio of the number of action potentials arriving at the terminal node A/9 to those generated in the initial node A/"o: 

transmission reliability coefficient: 1Z = (6) 

Ao 

The occurrence of a spike in the membrane potential Vi (t) is identified by an upward-crossing of the detection barrier 
at O.OmV ( Further details on the numerics can be found in Ref. [H|). The spike occurrences t\ with j = 1, ...,A/i 
where A/i indicates the number of spikes appeared on the specific node i define point processes Ui(t) = — 



A. Deterministic dynamics 

We start with the deterministic limit which is formally achieved in the limit A — > 00, using N = 10. In this 
case, the transmission reliability 1Z depends solely on the inter-nodal coupling strength k and exhibits distinguished 
transmission patterns [l5[ . 

The dependence of the transmission reliability 7Z on the coupling strength is depicted in FigJ5J For sub-threshold 
coupling, i.e. k < 0.067 mS/cm 2 the ratio equals zero and no spike propagation to the final node is observed. Contrary, 
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FIG. 3: (color online) Transmission reliability in presence of channel noise: The transmission reliability coefficient 1Z is plotted 
against the membranal size of the node of Ranvier for different strengths of the inter-nodal coupling (the coupling strengths 
are given in units of mS/cm 2 ) . The data confirms the result of noise-assisted spike propagation for sub-threshold coupling 
(see dashed red line) [HI]. Moreover, the minimum in 7L indicates the cross-over from a stochastic, uncorrelated spiking of the 
initial and final nodes towards a causal, noisy spike transmission from the initial to the final node. 
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FIG. 4: Spatial-temporal transmission patterns of spike propagation: the spatial-temporal evolution of the membrane dynamics 
at different axonal nodes of Ranvier is plotted for an intermediate coupling strength n = 0.08mS/cm 2 and various nodal 
membrane sizes: (a) A = 3 • I0 4 fim 2 (low channel noise level); (b) A = 100 iim 2 (intermediate channel noise level) and (c) 
A — 10 fim 2 (high channel noise level). Next to each panel, there is a color bar indicating the actual value of the membrane 
potential. The action potentials are created by a current stimulus at the initial node ("0"). In case of low noise level, cf. panel 
(a), the deterministic transmission pattern of 2 : 1 is clearly visible. In contrast, in the strong noise limit, due to the dominating 
spontaneous spiking, irregular transmission is found. 



for sufficiently large coupling parameter, i.e. k > 0.136 mS/cm 2 , each generated action potential propagates along the 
axon and arrives at the final node, i.e. TZ = 1. Discrete, rational transmission patterns k : I appear for intermediate 
values of the coupling parameter. 



B. Channel noise effects 



When considering finite sizes A of the nodes of Ranvier, the staircase-like dependence of the transmission reliability 
parameter TZ, which is depicted in Fig. [5] for the deterministic case (i.e. for the case of an infinite size of the nodes of 
Ranvier), turns in a continuous dependence (not shown). 

In the limit of strong intrinsic noise, i.e. for the finite size of the nodes of Ranvier approaching formally zero, the 
channel noise reigns the dynamics of the membrane potential of each node and the influence of the bi-linear coupling 
becomes negligible. Consequently, in this limit the causal relationship between spiking in the first and that one of 
the terminal node gets lost. The number of spike occurrences at the final node does not markedly depend on the 
actual value of the coupling parameter and the ratio TZ tends to a value independent of coupling strength k, cf. 
Fig. [3] Note, that in this case, the information transfer from an initial to a final node fades towards zero as the 
correlation between spikes in the initial and the final node diminishes. Moreover, the spontaneous spikes stimulates 
via the coupling the neighboring nodes and consequently a transmission to both sides occurs, cf. Fig. [4] (c). In the 
limit k — y 0, the initial and the terminal node are spiking independently. Consequently, 1Z is given by the ratio 
of the number of spontaneous spikes occurring within the non-stimulated, stochastic Hodgkin-Huxlcy dynamics and 
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the number of spikes occurring within the stochastic Hodgkin- Huxley dynamics complemented by a constant current 
stimulus of / = 12 /iA/cm 2 . Surely, the causal relationship is lost in this case. This becomes evident when considering 
the cross-correlation between the spiking of the initial and that of the final node (see below). 

With increasing size of the nodes of the Ranvier A, i.e. with decreasing channel noise strength, the transmission 
becomes more regular, cf. Fig. [4] (a) and (b). The causal relationship is restored in this limit and the transmission 
coefficient 1Z tends to its deterministic limit. This cross-over from the stochastic non-causal firing regime to the regime 
of noisy spike transmission depends on the coupling strength. For smaller coupling strengths, this cross-over occurs 
at larger nodal membrane sizes A, cf. Fig. [3J 

Starting out from the weak noise limit, an increase in the noise strength results in a noise-assisted spike propagation 
phenomenon for sub-threshold coupling as pointed out before with Ref. [l5j . The transmission reliability exhibits a 
maximum, indicating an optimal, noise assisted spike propagation, cf. Fig. [3] In this case the coupling between 
the nodes does not result in an efficient propagation of the spikes and the presence of intrinsic noise is necessary to 
overcome the threshold for excitation. However, for the supra-threshold coupling, the channel noise leads to noise- 
induced propagation failures and the transmission reliability coefficient 1Z firstly decreases with increasing noise level, 
cf. Fig. [3] The observed increase of 1Z for strong intrinsic noise is attributed to the cross-over to the stochastic, 
non-causal firing regime accompanied by uncorrelated spiking in the initial and final nodes. 

In order to analyze the spike correlation between the initial and final node, we firstly segmented the point processes 
uo(t) and ug(t + t) in segments of width At. For At smaller than the refractory time, there is either no spike or 
one spike observable in each segment. For our analysis, we chose At = 1.5 ms. Secondly, we determine the number 
of spike coincidences A/0,9 (t) between the initial node and the final node in the segments of width At, i.e. the spike 
coincidences between the two point processes Uo(t) and ug(t + r): 

A/o,9(T)=At / dtfo(t)f 9 (t + T), (7) 

Jo 

where the point process Ui(t) is approximated by 

/t+At 
dt'u t (t')/At for i = 0,9. (8) 

In Eq. [Jj T denotes the total integration time. Note, that the above definition corresponds to a cross-correlation 
measure. Thirdly, we relate this number A/0,9 (t) to the number of initiated spikes at node "0", i.e. Ao, and the 
bin-width At. We obtain the probability density Co,9(t): 

Note, that due to the periodic spike initiation at the first node, Co, 9(7") is periodic, cf. Fig. [SJ The normalization to 
the total number of initiated spikes Ao ensures that the integral of Co g(r) over one period results in the transmission 
coefficient 1Z defined by Eq. ©. 

A sharp peak in Co. 9(7") indicates a high correlation. For supra-threshold coupling and weak intrinsic noise, the 
correlation measure Co,g(r) exhibits a sharp peak for a r- value, which corresponds to the propagation time modulo 
the period of the spiking in the initial node, cf. Fig. [5] (b). With increasing noise level, the width of the peak 
increases and the height of the peak compared to the background level decreases until the total disappearance of the 
peak in the background showing the lost of any causal correlation. Interestingly, for the sub-threshold case, not only a 
broadening and flattening of the peak can be observed with increasing noise level, but also a shift of the peak towards 
smaller r- values showing a speed up of the transmission, cf. Fig. [S] (a) . The latter can be explained by the same line 
of reasoning as drawn by the effect of anticipated synchronization in Ref. [22| . 

IV. CONCLUSION 

We numerically studied the saltatory spike propagation along a myelinated axon within a mult i- compartment al, 
stochastic Hodgkin-Huxley modeling. The channel noise affecting the dynamics of the bi-linearly coupled nodes 
of Ranvier is originated in the random ion channel gating. Upon analyzing the spike propagation in terms of 
transmission reliability, i.e. the ratio of the number of spike observed in the terminal node to the number of spikes 
initiated in the first node, we found a reduction of the transmission reliability with increasing channel noise level 
for supra-threshold coupling. This is due to the noise induced propagation failures. A further increase of the 
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FIG. 5: (Color online) Initial-final-node spiking correlation: The spike correlation Cq,9(t) between the spiking at the initial 
node and that of the final node, cf. Eq. (J5J), is plotted for different coupling strengths k and membrane sizes A: In panel (a) 
for sub-threshold coupling k — 0.066 mS/cm 2 and membrane sizes A = 10 /zm 2 (blue solid line), A = 800 /im 2 (black dotted 
line), _4 = 3-10 3 /zm 2 ( green dash-dotted line) and A= 5 ■ 10 5 fim 2 (red dashed line); In panel (b) for supra-threshold coupling 
k — 0.08mS/cm 2 and membrane sizes A = 10 /im 2 (blue solid line), A = 100 /im 2 (black dashed line), .4 = 3- 10 4 /j,m 2 (red 
dash-dotted line). For weak intrinsic noise the sharp peak indicates a strong correlation. In contrast, in the strong noise limit 
the causal relationship is lost and the correlation function does not exhibit a peak. Interestingly, for the sub-threshold case 
a shift of the peak towards shorter times is observed, i.e. there is a speed up of the signal transmission with increasing noise 
level. 



channel noise level leads to the total loss of the spiking correlation between the first and the last node of Ranvier 
of the axonal chain. In case of sub-threshold coupling and not so high channel noise level, the channel noise can 
constructively contribute to the spike propagation and the effect of noise-assisted spike propagation is recovered. 
Both the transmission reliability as well as the propagation speed can increase with increasing channel noise level. 
This observed behavior is quite in spirit of the stochastic resonance phenomenon @-Q with an intrinsic noise source 
[lol. [Tl]| with the inter- nodal coupling playing the role of an nodal stimulus. 

The authors like to applaud and warmly thank Lutz Schimansky-Geier for his continuous engagement in furthering 
stochastic physics within the statistical physics community worldwide and for his many elucidative discussions with 
us in pursuing stochastic physics. He is still young and strong enough to appreciate and to contribute great science. 
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